Analysis of static and spherically-symmetric solutions in NDL theory of gravitation 
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We investigate a static solution with spherical symmetry of a recently proposed field theory of 
gravitation. In this so-called NDL theory, matter interacts with gravity in accordance with the 
Weak Equivalence Principle, while gravitons have a nonfinear self-interaction. It is shown that the 
predictions of NDL agree with those of General Relativity in the three classic tests. However, there 
are potential difi^erences in the strong- field limit, which we illustrate by proving that this theory 
does not allow the existence of static and spherically symmetric black holes. 



I. INTRODUCTION 

The Equivalence Principle played a fundamental role in 
the development of General Relativity (GR) and is at the 
heart of the idea that spacetime is curved However, a 
distinction must be made between the Weak Equivalence 
Principle (WEP) and the Strong Equivalence Principle 
(SEP). In a nutshell, WEP states that all matter fields 
(except the gravitational field) interact with gravity in 
the same way (the so-called "universal coupling"). SEP 
states, among other things, that all fields (including grav- 
ity) interact with gravity in the same way. Let us remark 
that while the other aspects of SEP (i.e. the absence of 
preferred-frame and preferred-location effects) have been 
tested in several instances Q , there is no conclusive evi- 
dence about the way that gravity couples to gravity. This 
is precisely the area explored by a recently presented field 
theory of gravitation, called NDL by the authors of 
This theory is a field-theoretical description of gravity (in 
the spirit of Feynman 1^ and Deser | ) in which SEP is 
violated from the beginning: gravity does not couple to 
itself in the same way it couples to other fields. 

Several features of this theory have been studied in 
lH, |[ 1^. To date, the most important prediction of NDL 
theory is that gravitational waves do not propagate at the 
same speed as electromagnetic waves a phenomenon 
related to the self-coupling of gravity in this theory. In 
fact, it can be shown that gravitons in NDL theory follow 
an effective metric, different to the background metric 
1^, that depends on the abovementioned self-coupling. 
Indeed, this is a feature of any nonlinear theory: the 
influence of the effective geometry on nonlinear photons 
has been extensively analyzed in ^. 

To go further with the analysis of the predictions of 
NDL theory, in this article we shall study some aspects of 
the static and spherically symmetric solution obtained in 
1^ . We begin by giving a short summary of NDL theory 
in Sectjl|. The spherically symmetric and static solution 



investigate the effective potential felt by particles moving 
in this solution. Armed with the effective potential, we 
compare the predictions of NDL with those of GR in the 
three classic tests. We also investigate if the spherically 
symmetric and static solution can describe a black hole 
analog to that of Schwarszchild in GR, and analyze the 
effective metric felt by gravitons. We close in Sect.^ with 
some comments regarding the nature of the singularities 
appearing in the solution and prospects for future work, 
get rid of eqn numbers 



II. SUMMARY OF NDL THEORY 

A detailed presentation of NDL theory was given in 
Q . Here we shall list the most important features of the 
theory, following Ref.[||: 

• The gravitational interaction is represented by a 
symmetric tensor ipf^^ that obeys a nonlinear equa- 
tion of motion. 

• Matter (but not gravity) couples to gravity through 
the metric g^i, — j^i, + ip^u, where 7^,^ is the flat 
background metric. 

• The self-interaction of gravity, given by a nonlin- 
ear Lagrangian, breaks the interpretation of grav- 
ity as a universal geometric phenomenon. That 
is, all particles (except gravitons) move follow- 
ing geodesic of g^j^. Gravitons instead move on 
geodesies of an effective m etric, the expression of 
which we shall give in Sect. dlVq ). 



is presented in Sect.HL We study in Sect, IV some aspects 



of gravitation around compact objects. Specifically, we 



It is convenient to define the tensor Fapfj_ (called the grav- 
itational field), in terms of (^^^ as follows [ |l^ : 

FaPt, = 2 (VM[a;/3] + F[all3]tj) , 

where the covariant derivative is constructed with the 
background metric 7^^. Indices are raised and lowered 
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with that metric also, and 

To construct a nonlinear theory of the gravitational field 
Ffii,a with the correct weak field limit, we assume that 
the interaction of gravity with itself is described by a 
functional oi A — B, where the scalars A and B are given 
by: 

A = Fa/3;.i^"^^, B = F^F^. 

From the action 

j d-^x y/^ C{A- B), 

(where 7 is the determinant of the background metric) 
and using that fact that Ca = —J~-b (where Cx is the 
derivative of the Lagrangian with respect to X), we get 
the equations of motion 



-7CaF 



0. 



(1) 



In the next section we shall study a particular solution 
of these equations of motions for a given choice of the 
Lagrangian. 



III. THE STATIC AND SPHERICALLY 
SYMMETRIC SOLUTION 

In what follows, we will restrict our study to a special 
Lagrangian studied previously in |Q , and inspired by the 
Born-Infeld Lagrangian : 



C{A -B) = — 

K 



1-^-1 
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(2) 



where k is Einstein's constant. The parameter b, with 
dimensions of length^^, is undetermined at this point. 
Notice that its value should be large enough for the series 
expansion of Eqn. (^) to be in agreement with the weak- 
field limit. 

We are interested in static and spherically symmetric 
solutions in a Minkowskian background. Consequently, 
the only nonzero components of the field ip^^i, are 



(pQQ = ^i{r), (fill 



and the corresponding nonzero components of the tensor 
FafSfj. are given by 



100 



r 



F122 = 2^^^ ~ 



^133 = ^122 sin^ 9. 
The only nonzero trace component, i^i, is given by 

Fl=^i' - -V. 

r 

From the equations of motion (|l]) and the expression of 
the Lagrangian we get only two nontrivial equations: 



fj.'r — ly — 

The first equation can be easily integrated, and the result 
is 



v{r) = e — 
r 



-1/2 



where we have defined 



C 



(3) 



(4) 



and C is an integration constant. From the second equa- 
tion. 



dr + const. 



Noting that the function ^{r) is defined only for r > tq, 
we can write using Eqn. (0) , 



- e C 



dr 



const. 



This integral can be written in terms of an elliptic integral 
of the first kind using the identity 



dx 



1 



valid for u > (3 > 0, and 



F{X,Y), 



X = arccos 

Consequently, /i(r) is given by 
e C 



Y 



V2 



ro 



y/a^ + /32 
i^(arccos(ro/r), 1/V2) -|- const. (5) 



In order to determine the value of the constant, we im- 
pose that spacetime be Minkowskian for large r, which 
means that //(r) — > in that limit. Then, 



\b\C 



f^(arccos(ro/r), 1/^2) - F(7r/2, 1/V2) 

(6) 



3 



The value of the constant C and the sign of e have 
not been determined up to now: they are dictated by 
the weak field limit of the solution, which must coincide 
with that of Schwarszchild. The radial component of the 
Schwarszchild metric (we are setting G = c— 1) is 



grr{r) 



-1 



In our case, grrir) 
limit we get 



2M 

r 

i/(r), and in the weak field 



9rr{r) 



C 



(7) 



and therefore conclude that e = +1 and C — 2M . 

Thus, the general expression for the static, spherically 
symmetric and asymptotically flat spacetime in NDL the- 
ory previously derived in |2| , is given by the metric 



ds^ = [1 + ^(r)] dt^ - [1 + v{r)] dr^ 



(8) 



with fi{r) and iy{r) given by Eqns. (^) and respec- 
tively. In the following sections we shall study some prop- 
erties of this solution. 



IV. GRAVITATIONAL PHYSICS AROUND 
COMPACT OBJECTS IN NDL THEORY 

Here the predictions of NDL theory are compared with 
those of GR. Using the static and spherically symmetric 
solution, we will consider first the three classical tests of 
GR in the framework of NDL theory. Then we shall see if 
NDL theory can describe black hole configurations analog 
to that of Schwarszchild. The following series expansions 
will be necessary in the ensuing sections. 



9uir) = l-- 
r 



1 



10 V r / 24 V r 



, (9) 



with Tq ~ 2M/\b\. The inverse of gtt and g^r are given 
by 

g^,\r) = 1 + ^ + ^ + ^+ 0(r-4), (10) 



where the dot indicates a derivative with respect to the 
proper time (or an affine parameter in the case of null 
geodesies), and 9 = tt/2. From the equation for the in- 
terval it follows that 



gut + grr-r + 



(13) 



where is 1 (0) for timelike (null) geodesies. Rearranging 
Eqns.(p^) and ( p^ we obtain 

f2 -I- V{r) = E^, (14) 

where the effective potential V{r) is given by 



V{r) 



1 



grr{r) 



1 



,gtt{r)grr{r) 



1 



with grr{T) and gttif) defined in Eqn.(g). L and E 
are the angular momentum and energy per unit mass 
respectively. This expression for V{r) reduces to the 
Schwarszchild case of GR, 



L2 



(15) 



upon substitution of the Schwarszchild metric. 



1. The perihelion shift of Mercury 

One of the triumphs of GR was its ability to correctly 
predict the perihelion of Mercury. Wc shall show here 
that NDL concurs with this prediction. The values 



3km, r+ w 7.0 x lO^km, r_ w 5.0 x lO'^km, 



will be used below, where (r_) is the aphelion (peri- 
helion) of Mercury. With these values, we see that 



|6P 



10 



-28 



This equation and the discussion following Eqn.(^ en- 
sure that we can keep only the first term in the expansion 
in Ecin.(ph. Notice also that 



;'M = -i + 7-7j + 7T + o(r-4)- (11) 



A. Motion of particles 

We now analyze the geodesies of photons and particles 
with nonzero mass in the geometry described by Eq.(^. 
The motion can be studied using the effective potential 

|l^ . For a static and spherically symmetric geometry 
there are two constants of motion along a geodesic, which 
we designate by E and L: 



gut = E, 



(12) 



r± 



10 



-29 



and consequently from the expansions in Eqns. (10) and 
( |rT| ) the fourth and higher orders terms can be neglected. 
The effective potential in the case under consideration 
can thus be approximated by 



V{r) = k\l 



I? 



+ '^{k~E^). 



(16) 

We see that this approximate expression contains new 
terms not present in the exact expression given by 
Eqn. (p^) , even for the massless case. These new terms 
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signal potential differences between NDL theory and GR 
in the weak field regime. Let us try to estimate the mag- 
nitude of these terms for the case of Mercury, in which 
C^/L^ « 10"'', so the term involving this factor can be 
safely neglected. We now need an estimate for which 
may be obtained using Eqn.(|lJ) evaluated at r+ and r_ 
(the values of the radius where r = 0). The result can 
be written as a system of two equations in the unknowns 
and L^. In particular, 

^2 ^ 9tt{r+)gu{r-) {rl - rl) 



r+9tt{r-) - r±guir+) 

It follows that E^ — 1 K 10^'' and so the last term in 
Eqn.(|l^) is negligible for fc = 1. Consequently the ef- 
fective potential for NDL theory coincides with that of 
Schwarszchild. Therefore, NDL agrees with GR in the 
prediction for the perihelion of Mercury. 



2. The deflection of light by the Sun 

To derive the equation that governs the deflection of 
light rays by the Sun, we use the fact that 



Rq 



4.23 X 10" 



(17) 



Combining this with the discussion following Eqn.(|^), it 
follows that the effective potential for the deflection of 
light is given by Eqn.([l6|) with fc = 0. That is, 



V{r) 



L2 

,2 



1 



E' 



\ r / 

Consequently, the equation of motion takes the form 



C 



^ 1-- = 



1 + I7 



E. 



From Eqn.(|T^) we have that (C/r)^ < 10^^^, which is 
negligible. Thus it is clear that the equation of motion 
for photons in NDL coincides with that of Schwarszchild. 



3. Time delay of light 

Let us sketch the derivation of this effect in NDL the- 
ory, adopting the usual procedure for GR (after [p^). 
Consider the path of a light ray with 9 = tt/2 in the 
metric given by Eqn.(^): 



gttdf + grrdr'^ 
This equation can be rewritten as 

2 



dr = dr^ 



dr 



9tt 9rr9tt 



Using the expansions given in Eqns.(|^) and (|o|) we get 

dt^ = dr"^ 



dr 



2C 
r 



which coincides with the equation for the GR case (see 
JlSf , pp. 204). Consequently, also in this case NDL the- 
ory prediction agrees with that of GR. 



B. Black holes in NDL 

We now investigate the existence of horizons in the 
metric given by Eqn.(||). The position of the putative 
horizons is given by the condition gtt{rh) = 0. This im- 
plies 



i^(arccos(ro/r;,), 1/V2) - F{Tr/2, 1/V2) 



-x/2^. 
C 



By standard manipulations (see for instance p4||) we ob- 
tain 



rh = ro 



F(7r/2,1/V2) 



C\b\ 



(18) 



where the function cn(a;) is the cosine-amplitude. This 
expression provieds the value of in terms of C = 2M 
and the constant |6|. Note that rh > tq always. Con- 
sequently, the region where the geometry is not defined 
(r < ro) is always inside the surface r = rh- 

To determine whether Eqn.(|l8|) defines a horizon or a 
singular surface, we can compute the components of the 
Riemann tensor calculated in the tetrad system 



ujg = r, ^0 = rsmt'. 

Here we will give the expression for only one of the com- 
ponents, which will be enough to illustrate the typical 
behaviour: 



Rt$t0 = 



M'(r) 



2r{l + v{r)){l + ^i{r))' 



(19) 



Clearly this component of the Riemann tensor in the 
tetrad system depends on the product {gtt{r)grr{'r))~^ 
which diverges at r = r;i because <7tt(f) is null at that 
point. Consequently, we would have a naked singular 
surface instead of a horizon We could attempt to 

cure this divergence by imposing that grrif) diverges at 
the horizon, in the hope of getting a finite and nonzero 
value for the denominator of Eqn.(|l^). Note however 
that grr{f) only diverges at r = tq, so it would be nec- 
essary to impose that rh = ro to produce a horizon. To 
examine the behaviour of the product for rh = ro, we 
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expand the metric functions near r = tq. It is convenient 
to change variables to r = tq + R. In terms of this new 
variable, the significant terms in the series expansions are 



C 



(ro + R) 



9C / [R 

arctan W — 

8ro \y ro 




gttiro+R) ~ 1 



The only term that leads to a divergence in the limit 
i? — > is the term proportional to 1/ appearing in the 
second expression. Taking the product gttiTo+R)grr{ro + 
R) in the limit R ^ (and using that arctan(a;) « a; — ^ 
for small x), it is clear that the divergence cannot be 
eliminated. The analysis of the geometry seen by matter 
then implies that static and spherically symmetric black 
holes cannot be described by NDL theory: the would-be 
horizon r = rh is in fact a naked singularity. 



C. Effective metric for gravitons 

We have shown that there are no black holes for mat- 
ter in NDL theory. However, there may be a black hole 
configuration for gravitons. As previously mentioned, in 
NDL theory gravitons interact nonlinearly with them- 
selves. Consequently, the path of these particles is not 
governed by the background metric but by an effective 
metric given by 



Pfiu = It^u + A^ 



with 



An. = 2 



c 



AA 



Ca 



and £{A—B) is an arbitrary Lagrangian. A short calcula- 
tion using the Born-Infeld-like Lagrangian from Eqn.(|^) 
and the metric given by Eqn.(0) shows that 



Pooir) 



1 



v{rf 



with pii{r) — —poo{r), while p22(f) and pssir) are as 
in Minkowski spacetime. The result for the compo- 
nent of the Riemann tensor given by Eqn.(^9|) calcu- 
lated with this metric indicates that there is a (naked) 
singularity seen only by gravitons at r = rp. Subse- 
quently, the effective metric for gravitons cannot describe 
a Schwarszchild black hole. 



Note however that there are potential disagreements in 
situations in which gravity is strong. For instance, for 
compact objects like neutron stars, for which typically 
M / Rk, 1Q^Mq/Rq, the convergence of the series expan- 
sions Eqns.(|9|)-(lTl|) will be much slower, and observable 
effects would appear. Another instance in which the two 
theories differ in their predictions is in the existence of 
black holes. As we have shown, there are no solutions in 
NDL theory that describe static and spherically symmet- 
ric black holes. In fact, from the analysis in Sect. (Ill), we 
conclude that the two possible cases in the metric seen 
by matter lead to naked singularities. 

Let us remark that there are two different types of sin- 
gularities appearing in the problem, which we will call 
geometrical and physical. For rg r^, we encountered 
two singularities: one located at r — rh and another 
located at r — rg. We can call the first singularity geo- 
metrical because the quantities related to the geometry 
(like the effective potential and the Riemann tensor in 
tetrads) diverge at r = r/i, but the field is finite 

at that point (see Eqn.(II])). The energy density of the 
field is also finite: its energy-momentum tensor is given 
by @ 

T^, = -C{ABh^. + 2[2F;'^F,^p + Fo,f3^F''t-F^F, 

From this expression, and using the equations of motion, 
we get for the energy density, 

2 / 

^0" = -13^ + 2^(/i - ly'r) + 2nj, 

which is finite at r = r^. On the other hand, at r = rg 
the quantities related to the geometry do not diverge, 
but there is a divergence of the field and of the energy 
density. 

The divergences of the geometrical and field quantities 
are not compatible for rg = rh either. We expect that 
the divergences of the quantities related to the field to 
be physical. The fact that the geometry does not display 
the same pattern of divergences as that of the field sug- 
gests that some modification should be introduced to the 
theory. A possible alternative would be to assume that 
the (unobservable) background is not flat spacetime but a 
curved geometry (for instance, de Sitter spacetime). An- 
other assumption that must be checked is the postulated 
coupling between matter and gravity (given by the met- 
ric g^i, — "ffj_,y + i^^u)- Work in this direction is currently 
under way. 



V. DISCUSSION 

Our analysis revealed that NDL agrees with GR for 
the three classical tests performed in the weak field limit. 
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